function I = Kc(Itot)
IR = Itot(:, :, 1);
IG = Itot(:, :, 2);
IB = Itot(:, :, 3);

[M, N] = size(IR);
I(1 : 4 : M, 1 : 4 : N) = IR(1 : 4 : M, 1 : 4 : N);
I(1 : 4 : M, 2 : 4 : N) = IR(1 : 4 : M, 2 : 4 : N);
I(2 : 4 : M, 1 : 4 : N) = IR(2 : 4 : M, 1 : 4 : N);
I(2 : 4 : M, 2 : 4 : N) = IR(2 : 4 : M, 2 : 4 : N);


I(1 : 4 : M, 3 : 4 : N) = IG(1 : 4 : M, 3 : 4 : N);
I(1 : 4 : M, 4 : 4 : N) = IG(1 : 4 : M, 4 : 4 : N);
I(2 : 4 : M, 3 : 4 : N) = IG(2 : 4 : M, 3 : 4 : N);
I(2 : 4 : M, 4 : 4 : N) = IG(2 : 4 : M, 4 : 4 : N);


I(3 : 4 : M, 1 : 4 : N) = IG(3 : 4 : M, 1 : 4 : N);
I(3 : 4 : M, 2 : 4 : N) = IG(3 : 4 : M, 2 : 4 : N);
I(4 : 4 : M, 1 : 4 : N) = IG(4 : 4 : M, 1 : 4 : N);
I(4 : 4 : M, 2 : 4 : N) = IG(4 : 4 : M, 2 : 4 : N);


I(3 : 4 : M, 3 : 4 : N) = IB(3 : 4 : M, 3 : 4 : N);
I(3 : 4 : M, 4 : 4 : N) = IB(3 : 4 : M, 4 : 4 : N);
I(4 : 4 : M, 3 : 4 : N) = IB(4 : 4 : M, 3 : 4 : N);
I(4 : 4 : M, 4 : 4 : N) = IB(4 : 4 : M, 4 : 4 : N);


end